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ABSTRACT 

We present the Spitzer Space Telescope Infrared Array Camera (IRAC) observa- 
tions for a sample of local elliptical galaxies to study later stages of AGN activities. 
A sample of 36 elliptical galaxies is selected from the Palomar spectroscopic survey. 
We detect nuclear non-stellar infrared emission in 9 of them. There is unambiguous 
evidence of circumnuclear dust in these 9 galaxies in their optical images. We also 
find a remarkable correlation between the infrared excess emission and the nuclear 
radio/X-ray emission, suggesting that infrared excess emission is tightly related with 
nuclear activity. Possible origin of infrared excess emission from hot dust heated by 
the central AGN is supported by spectral indices of IR excess emission. 



Key words: galaxies: active 
- infrared: galaxies 



galaxies: elliptical and lenticular, cD - galaxies: nuclei 



O 1 1 INTRODUCTION 



The tight correlation between the velocity dispersion of the 
bulge of the host galaxy and the black hole mass suggests 
that supermassive black hole (SMBH) is an ubiquitous com- 
ponent in elliptical galaxies (Gebhardt et al. 2000; Ferrarese 
k, Merritt 2000; Kormendy & Richstone, 1995; Magorrian et 
al. 1998). Activities of SMBHs are well known for elliptical 
galaxies at high redshifts, these systems are host galaxies for 
most radio-loud and the bright quasars (Hutchings & Mor- 
ris 1995; Bahcall, Kirkakos & Schneider 1996; Falomo et al. 
2005). However, most SMBHs in nearby ellipticals are no 
longer very active. Ho et al. (1995, 1997) found that, in the 
Palomar spectroscopic survey of nearby galaxies, only about 
50% of ellipticals show detectable emission-line nuclei, most 
of which belong to be low-ionization nuclear emission-line 
regions (LINERs; Heckman 1980). 

Mulit-wavelength band surveys in elliptical galaxies 
have been carried out to search for nuclear emission from 
LINERs and other low-luminosity active galactic nuclei 
(LLAGNs) to study the true ionization mechanism. Such ob- 
servations were performed in optical (Ho, et al. 1995, 1997); 
X-ray (Terashima, et al. 2000, 2002; Flohic et al. 2006), ra- 
dio (Nagar et al. 2002, 2005; Filho et al. 2006), and ultravi- 
olet band (Maoz, et al. 1995; 1996). None of these observa- 
tion were conclusive in determining the nature of these ob- 
jects. A widely accepted perspective is that LLAGNs display 
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remarkably different properties with their high-luminosity 
counterparts such as quasars and Seyferts. The optical-to- 
ultraviolet "big blue bump", a typical Spectral Energy Dis- 
tribution (SED) of high-luminosity AGNs, is weak or absent 
in LLAGNs (Ho, 1999) . At the same time, an anticorrelation 
between radio-loudness and Eddington ratio for AGNs is re- 
ported and confirmed by recent studies (Ho, 2002; Greene 
et al. 2006; Sikora et al. 2007). In the infrared band, the 
composite SEDs of LLAGNs display a mid-infrared peak 
(Ho, 2008). Ho (2008) concludes that LLAGNs have an ra- 
diatively inefficient accretion flow (RIAF) with a truncated 
thin disk. Maoz (2007), however, found that SEDs for a sam- 
ple of relatively unobscured low luminosity LINERs show 
no difference with that for higher luminosity AGNs. A sim- 
ilar conclusion was obtained recently by Dudik et al. (2009) 
that [NeV]24/im /[OIV]26/im mid-infrared line flux ratio for 
LLANGs similar with standard AGN, arguing against a UV- 
to-optical deficiency due to inefficient accretion in LLANGs. 

On the other hand, recent observations show that the 
circumnuclear region of elliptical galaxies is also compli- 
cated. Elliptical galaxies are thought to contain only old stel- 
lar population and hot gas. This picture has been challenged 
by recent observations. Shields (1991) detected warm gas 
(T ~ W 4 K) in many elliptical galaxies (see also Macchetto 
et al. 1996). Recent neutral hydrogen observations reveal 
substantial amount of neutral hydrogen gas in many early- 
type galaxies (Morganti et al. 2006; Noordermeer 2006). 
Moreover, even cool interstellar medium (ISM) including 
dust and molecular gas were detected in elliptical galax- 
ies (Knapp et al. 1985; Knapp et al. 1989; van Dokkum et 
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al. 1995; Wiklind et al.1995; Temi et al. 2004; Lauer et al. 
2005; Sage et al. 2007; Kaneda et al. 2008). The circumnu- 
clear cold interstellar media are mostly detected in elliptical 
galaxies with nuclear activities (Tran et al. 2001; Krajnovic 
& Jaffe, 2002; Xilouris & Papadakis, 2002; Simoes Lopes et 
al. 2007; Zhang et al. 2008). It is more likely that active 
elliptical galaxies possess circumnuclear dust feature. This 
fact challenges the RIAF accretion scenario, which assumes 
LLAGNs are not gas-starving, but have a low radiative effi- 
ciency. 

The Spitzer Space Telescope (Werner et al. 2004) with 
2" spatial resolution (Fazio et al. 2004) in mid-infrared 
(3.6 — 8.0fj,m) offers us a new approach to study both 
LLAGNs and their environment in ellipticals. Although con- 
tributions of photospheric emission from evolved stellar pop- 
ulation and hot dust in the circum-stellar envelopes of AGB 
stars form a considerable mid-infrared background in ellip- 
ticals (Athey & Bregman 2002; Xilouris et al. 2004; Temi et 
al. 2008), smooth spatial distribution of surface brightness 
of elliptical galaxies permit to separate non-stellar nuclear 
source from galaxy component. There were several detec- 
tions of mid-infrared core feature in previous studies (Pahre 
et al. 2004; Gu et al. 2007). This excess emission from AGN 
component can not be easily extracted from host galaxy ex- 
tended component due to poor spatial resolution. 

In this paper, we present results of MIR observation of 
a sample of nearby elliptical galaxies at mid-infrared wave- 
length with Infrared Array Camera (IRAC) and Multiband 
Imaging Photometer (MIPS) aboard the Spitzer Space Tele- 
scope in search for mid-infrared core, and studies of correla- 
tion between infrared core feature and other nuclear proper- 
ties of elliptical galaxies. This MIR observation permits to 
investigate the origin of nuclear infrared emission in ellipti- 
cals. 

This paper is organized as follows: the sample selection 
and data reduction are described in Section 2; the basic re- 
sults of the observation are presented in Section 3; the results 
are discussed in Section 4; the conclusion is given in Section 
5. 



2 SAMPLE SELECTION AND DATA 
REDUCTION 

A sample of 36 elliptical galaxies were selected from the 
Palomar optical spectroscopic survey for this study. The 
Palomar optical spectroscopy survey comprises all nearby 
galaxies brighter than Bt = 12.5mag in the northern hemi- 
sphere (Ho et al. 1995, 1997). This sample is statistical 
complete, and contains both galaxies with and without nu- 
clear activity while the latter ones serve as a control sample. 
Multi-band observations of nuclear region has been carried 
out and studied for a large fraction of this sample, making it 
suitable for our study of the infrared properties. Our sample 
includes 2 Seyferts, 15 LINERs, 5 transition objects and 15 
inactive galaxies. 

The IRAC Basic Calibrated Data (BCD) and MIPS 
Post Basic Calibrated Data (Post BCD) of these galaxies 
are downloaded from the archive of Spitzer Science Center. 
The IRAC BCD images were performed with basic image 
processing, including dark subtraction, detector lineariza- 
tion corrections, flat-field corrections, and flux calibration. 



Table 1. Global and Nuclear Properties. 



Galaxy 


D 


M B 




Nuclear 




Name 


(Mpc) 


(mag) 


Class 


Dust 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC315 


65.8 


-22.22 


LI. 9 


Disk/Ring 


3 


NGC410 


70.6 


-22.01 


T 






NGC777 


66.5 


-21.94 


S 






NGC821 


23.2 


-20.11 


A 


No Dust 


2 


NGC1052 


17.8 


-19.90 


LI. 9 


Lane/Chaotic 


5 


NGC2832 


91.6 


-22.24 


L2: 


No Dust 


1 


NGC3226 


23.4 


-19.40 


LI. 9 


Disk/Ring 


3 


NGC3377 


8.1 


-18.47 


A 


Lane/Chaotic 


1 


NGC3379 


8.1 


-19.36 


L2/T2: 


Disk/Ring 


2 


NGC3608 


23.4 


-20.16 


L2/S2: 


No Dust 


2 


NGC3610 


29.2 


-20.79 


A 


No Dust 


2 


NGC3640 


24.2 


-20.73 


A 


No Dust 


1 


NGC4125 


24.2 


-21.25 


T 


Disk/ring 


1 


NGC4168 


16.8 


-19.07 


S 


Lane/Chaotic 


1 


NGC4261 


35.1 


-21.37 


L2 


Disk/Ring 


3 


NGC4291 


29.4 


-20.09 


A 


No Dust 


2 


NGC4278 


9.7 


-18.96 


LI. 9 


Lane / Chaotic 


2 


NGC4374 


16.8 


-21.12 


L2 


Lane / Chaotic 


5 


NGC4406 


16.8 


-21.39 


A 


No Dust 


2 


NGC4473 


16.8 


-20.10 


A 


No Dust 


2 


NGC4552 


16.8 


-20.56 


T 


Lane / Chaotic 


2 


NGC4564 


16.8 


-19.17 


A 


No Dust 


2 


NGC4621 


16.8 


-20.60 


A 


No Dust 


1 


NGC4636 


17.0 


-20.72 


LI. 9 


Lane/Chaotic 


7 


NGC4649 


16.8 


-21.43 


A 


No Dust 


2 


NGC4660 


16.8 


-19.06 


A 


No Dust 


2 


NGC5077 


40.6 


-20.83 


LI. 9 


Lane/Chaotic 


6 


NGC5322 


31.6 


-21.46 


L2: 


Disk/Ring 


2 


NGC5557 


42.6 


-21.17 


A 


No Dust 


1 


NGC5576 


26.4 


-20.43 


A 


No Dust 


2 


NGC5813 


28.5 


-20.85 


L2 


Disk/Ring 


1 


NGC5831 


28.5 


-19.96 


A 


No Dust 


1 


NGC5846 


28.5 


-21.36 


T 


Lane/Chaotic 


1 


NGC5982 


38.7 


-20.89 


L2 


No Dust 


1 


NGC6482 


52.3 


-21.75 


T 






NGC7626 


45.6 


-21.23 


L2 


Dust/Ring 


3 



Notes — Col(4):Nuclear spectral type from the Palomar survey 
(Ho et al. 1997):L=LINER; S=Seyfert; T=Transition object; 
A=absorption-line nuclei (inactive). Col(5):Morphology of opti- 
cal circumnuclear dust, obtained from references listed in Col(6), 
including three types: nuclear dust disk or dust ring; dust lane or 
disorganized dust patch; no dust. Col(6):References — (1) Tran et 
al. 2001; (2) Lauer et al. 2005; (3) Gonzalez Delgado et al. 2008; 
(4) Lauer et al. 2007; (5) Ravindranath et al. 2001; (6) Rest et 
al. 2001; (7) van Dokkum & Franx, 1995. 

We further use the custom IDL software (Huang et al. 2004) 
to make the final mosaic image for each object. The absolute 
flux calibration for IRAC flux densities is better than 10% 
(Fazio et al. 2004). We adopt the AB magnitude system for 
magnitudes and colors throughout this paper. 

To obtain the MIR color distribution of each galaxy, we 
first cross- convolve each image by using the corresponding 
PSF0(Gordon et al. 2008, Tom Jarrett, private communi- 



1 http: / / dirty, as. arizona. edu /kgordon /mips/ conv_psfs / conv_psfs. html 
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cation). For example, the color difference between image at 
3.6/im and at 8.0/im was obtained by following steps, firstly: 

Image(3.6') = Image(3.6) ® PSF(8.0) (1) 

Image(8.0') = Image(8.0) ® PSF(3.6) (2) 

Where "®" means convolution, then: 

Color(3. 6-8.0) = (Image3.6' x A 3 .e) - (Image8.0' x A 8 .o)(3) 

Where A3.6 and Ag.o are PSF aperture correction factors for 
an infinite aperture. 

We used ellipse program in the ISOPHOT package of 
IRAF to perform the surface photometry of each galaxy. 
Hot pixels and foreground stars were identified by eye and 
masked out before isophotal fitting. 

The isophotal parameters, such as ellipticity and posi- 
tion angle, were measured at 3.6/im where the Signal-to- 
Noise (S/N) ratio is the highest. These parameters were 
then applied for the surface photometry at 4.5, 5.8 and 
8.0/ira. Considering the scattered light, we employed ex- 
tended sources aperture correction for calibration provided 
by Tom Jarret10. Since we particularly focus on the nuclear 
region, we also obtained the nuclear flux density from a cir- 
cular region within an aperture of 10" (for NGC 1052 and 
NGC 3226, with an aperture of 15") to extract the non- 
stellar excess emission. The size of such a central region is 
determined by the radial color distribution shown in Figure 
1. 



3 RESULTS 

Figure 1 shows color distributions of [3.6] — [4.5], [3.6] — [5.8], 
[3.6] — [8.0] for all elliptical galaxies in our sample. In 
10" < R < 40", the infrared colors for almost all galax- 
ies do not change significantly, and are generally consistent 
with photospheric emission of late-type stars with a minor 
contribution of hot dust in circum-stellar envelopes of AGB 
stars (Pahre et al. 2004; del Burgo, Carter & Sikkema, 2008; 
Temi, Brenghti & Mathews, 2008). However, the colors in 
the central region with R < 10" tell a different story: 9 out 
of 36 galaxies exhibit much redder colors. The [3.6] — [8.0] 
shows most color excess in the central region. Thus we use 
[3.6] — [8.0] colors to distinguish galaxies with significant ex- 
cess emission from normal galaxies. For galaxies with no 
significant color excess, the deviation of [3.6] — [8.0] color in 
central 10" region is only —1.44 ± 0.02. The remaining 9 
galaxies all have [3.6] - [8.0] color redder than -1.34, indicat- 
ing a redder central color above 3cr level. 

Figure 2 is the color-color diagram for galaxies with no 
color excess in the center. The dispersions of their IRAC 
colors are very small, with [3.6] — [4.5] = 0.61 ± 0.02, 
[3.6] - [5.8] = -1.11 ± 0.02, [3.6] - [8.0] = -1.44 ± 0.02, in 
agreement with colors of M-type star ( Pahre et al. 2004). 
This result supports that normal ellipticals are dominated 
by old stellar population. Moreover, LLAGNs and inactive 
galaxies show no systematic difference in mid-infrared color, 
which agrees with the nuclear stellar population analysis of 
LLAGNs in previous studies (Boisson et al. 2000; Ho et al. 

2 http: / /spider. ipac.caltech.edu/stafT/jarrett/irac/calibration/ 



Table 2. Multi-wavelength Nuclear Emissions. 
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cllllC 


(erg/s) 


[erg/s/ Hz) 


(erg/s) 


(erg/ s/ Hz) 


(erg/s/ H 


K 1 ) 








\PJ 


\yi 


NGC315 


39.38 


30.39 


41.64 


30.22 


30.27 


NGC410 


< 39.32 


< 27.78 










NGC777 


< 39.11 


< 27.90 










NGC821 











22.45 


NGC1052 


39.43 


29.14 


40.78 


29.53 


29.72 


NGC2832 


< 39.05 


< 28.18 








NGC3226 


38.93 


27.55 


39.62 






NGC3377 








28.04 


28.39 


NGC3379 


37.73 


< 25.89 


37.89 








NGC3608 


37.80 


< 26.99 


38.64 






NGC3610 














NGC3640 














NGC4125 


38.85 


< 26.85 


38.93 


29.70 


30.02 


NGC4168 


37.91 


27.01 







29.30 


NGC4261 


39.70 


29.65 


40.65 


29.07 


29.28 


NGC4291 














NGC4278 


38.88 


28.10 


39.96 


28.83 


29.27 


NGC4374 


39.03 


28.79 


39.58 


29.24 


29.54 


NGC4406 








28.57 


29.00 


NGC4473 














NGC4552 


38.05 


28.30 


39.41 


28.73 


29.20 


NGC4564 














NGC4621 














NGC4636 














NGC4649 








29.43 


29.52 


NGC4660 














NGC5077 


39.52 










NGC5322 


38.54 


28.18 




29.71 


30.03 


NGC5557 














NGC5576 








28.88 


29.20 


NGC5813 


38.35 


27.37 










NGC5831 












NGC5846 


38.32 


27.79 


39.54 








NGC5982 


38.55 









29.77 


NGC6482 




< 27.52 


39.40 








NGC7626 


38.52 


29.00 











Notes — Col(2): Luminosity of [O/iT] A5007 taken from Ho et al. 
(1997), except NGC 1052, NGC 4125 and NGC 5813, which are 
based on observations under non-photometric conditions. Ho et 
al. (2003b) give their updated Ha luminosities and we use them 
to derive their [OI 1 1]\5007 luminosities, assuming same ratios of 
[OIII]\5007/Ha from Ho et al.(1997). Col(3): Nuclear Luminos- 
ity at 15GHz, taken from Nagar et al.(2005), except NGC1052, 
which is taken from Kellermann et al.(1998). Col(4): Nuclear Lu- 
minosity at 2-10keV, taken from Gonzalez-Martin et al. (2006), 
except NGC 4278, which is taken from Terashima & Wilson 
(2003). Col(5) - (6): Total luminosities at 60(j,m and lOOfim from 
IRAS observation (Knapp et al. 1989) 

2003b). Zhang et al.(2008) studied the nuclear stellar pop- 
ulation for a sample of early-type galaxies that is highly 
overlapped with our sample, and they found no difference 
in stellar age distribution between LLAGNs and inactive 
galaxies in their sample. 

The mid-infrared emission from elliptical galaxies con- 
sists of two components: stellar and non-stellar emission 
(dust, AGN, etc.). The 3.6/xm emission is dominated by 
later-type stellar photospheric emission (Pahre et al. 2004; 
Temi et al. 2008), thus traces stellar mass distribution. The 
non-stellar component becomes significant at longer wave- 
length, which is shown clearly by the [3.6] — [8.0] color. To 
derive the flux density of excess non-stellar emission in in- 
frared core galaxies, we use the mean colors ([3.6]- [5.8], [3.6]- 
[8.0]) for central regions of non-core galaxies as a template 



4 Yuping Tang, Q.-S. Gu, J.-S. Huang, Y.-P. Wang 



j i j { i i i i i i i i I I I I I I I i M 



;ri -0.5 -- 

S-'-Ori | | j j 

I 



; i • ; • i i i i i i i ! i I ! i 



■ liliiii 



i i i ! j I j i 1 1 



20 

SMA(orcsec) 



Figure 1. Color distributions of all non-core galaxies (black) in 
comparison with core (red) galaxies, the x-axis refers to distances 
from the center along major axis, in units of arcsec. Core galaxies 
have obvious redder color in central several arcseconds, especially 
seen from 3.6 — 8.0fim. 



of the old stellar component. For five infrared core galaxies 
with visible excess emissions (NGC 315, NGC 1052, NGC 
3226, NGC 4261 and NGC 5322) at both 5.8/im and 8.0^m, 
by assuming that the shape of SED for the central excess 
emission as a power-law function, we are able to disentangle 
excess component from stellar emission by solving following 
equations: 



fexs,i ~T" fs. 
fstar 7 i I fstc 



- ftot 



fexs,i/fexs,j — (^/^j) 



(4) 
(5) 
(6) 

where subscript numbers i/j = 1,2,3 correspond to differ- 
ent wavelength Ai,2,3 = 3.6,5.8,8.0, i^/j are corresponding 
frequency of Xi/j, f e xs,i/j an( i fstar,i/j are nu x densities of 
core emission and stellar emission at different bands, respec- 
tively. f eX s,i/j-> fstar,i/j an d a are set as variables, with a to- 
tal number of 7, equal to the number of equations. Ri,j are 
flux ratios of emission for old stellar population, obtained 
by averaging central colors of all non-core galaxies, ftot.i are 
flux densities of stellar emission at different bands, extracted 
from central region within an aperture of 15" for NGC 1052/ 
NGC 3226 and 10" for other core galaxies. The size of aper- 
ture is determined by two factors, it should be large enough 
to include the whole infrared core structure, and as small 
as possible to reduce the influence from offset of zero-point 
caused by uncertainty in true nuclear stellar color, this influ- 
ence bring larger uncertainty at shorter wavelength, where 
dilution by stellar emission is more serious. The special res- 
olution of IRAC is 2". In [3.6 - 8.0] distribution shown in 
Figure 1, except NGC 1052, NGC 3226 and NGC 4278, the 
sizes of red core structures for all the rest six galaxies are 
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Figure 2. Central colors of non-core galaxies within 10", stars 
are LLAGNs while triangles are inactive galaxies. 



smaller or close to 10 arcsec. Furthermore, due to the con- 
volution procedure mentioned in Section 2, the sizes of core 
structure in Figure 1 appear larger than their real values. 
Therefore, we consider 10 arcsec a proper size for the aper- 
ture to extract nuclear excess emission. For NGC 1052 and 
NGC 3226, 15" aperture is adopted, NGC4278 show a red- 
der color throughout the galaxy and a vague core structure, 
we simply use 10" aperture to derive a "nuclear" flux in this 
object. For such apertures, the propagated relative errors 
of excess emission due to uncertainties of Rij are several to 
about ten percent. 

Flux densities and luminosities obtained through above 
approachs are listed in Table 3 and 4. In Table 3, the errors 
of total IRAC and MIPS flux densities of whole galaxies are 
simply set as 10% (Fazio et al. 2004; Rieke et al. 2004), which 
is an estimation for the uncertainties of absolute calibration, 
in contrast, uncertainties derived from image statistics are 
neglibile. The uncerntainties of excess emission and spec- 
tral indices are estimated from deviation of excess emissions 
extracted from ±5" apertures (10"/15"/20" for NGC 1052 
and NGC 3226, 5"/10"/15" for other galaxies). The propor- 
tion of excess emission at 3.6 is very small, contributes less 
than 5% even in the strong infrared core galaxy, NGC 1052, 
and less than 0.5% in NGC 4261. Hence it is reasonable to 
assume all 3.6/im flux densities in the other four galaxies 
(NGC 4125, NGC 4278, NGC 4374 and NGC 5077) are pro- 
duced by stellar component, with this assumption we obtain 
the excess flux densities at 8.0/im for these four galaxies, 
since their excess emissions at short wavelength, even if ex- 
ist, will be too weak to be detected. The excess emissions 
are generally small comparing with stellar emission except 
in NGC 1052, where 8.0/im excess emission contributes to 
nearly 50 percents of total emission. 

Figure 3 shows 8. 0pm residual images of nine infrared 
core galaxies, which are obtained by utilizing 3.6/im image 
and flux ratio /Z3.6.8.0 appearing in equation (5) to remove 
the contribution of underlying stellar population. The 8/im 
excess emissions in four galaxies (NGC 315, NGC 1052, NGC 
4261 and NGC 5322) show point-like structure with ring- 
like feature of 8.0/im Point Spread Function (PSF), and the 
other five galaxies show extended excess emission, indicat- 
ing off-nuclear sources of 8fim excess. NGC 1052 show a 
substructure on the right of the center. This is due to the 
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Figure 3. Residual images of infrared core galaxies at 8fim in central 24", x24" regions, with the background from old stellar population 
removed. Point-like structure can be seen in four galaxies: NGC315, NGC1052, NGC4261 and NGC5322, while the rest galaxies show 
extended excess emission at 8fim. The substructure in NGC 1052 is a bandwidth effect. 



bandwidth effect and was masked for measuring photome- 
try. 

To show the contribution of extended emission, in Fig- 
ure 4 we compare nuclear 8/j,m excess luminosity with 
the compactness of excess emission, quantified by L(< 
5")/ Lex, total, the proportion of nuclear excess luminosi- 
ties within 5" aperture to total extended luminosities. For 
point-like sources, the surface brightness profile within 5" 
are uniform and consistent with 8/im PSF. The propor- 
tion of central emission within 5" are generally higher than 
50% and are higher than 80% for point-like sources. The 
compactness decreases as nuclear luminosity decreases. Ex- 
tended emissions become considerable only for sources with 
nuclear luminosities lower than 3 x 10 27 erg / s/ Hz. 



4 DISCUSSION 

The excess non-stellar emission may originate from the cen- 
tral AGNs or nuclear hot dust heated by AGN. PAHs emis- 
sion feature could also contribute to excess emission at 
8.0 \im, this is exactly the case for NGC 1052, where 7.7fim 
PAHs emission feature has been detected (Kaneda et al. 

3 IRAC Data Handbook: http://ssc.spitzer.caltech.edu/irac/dh/ 



2008), though this feature by itself is not able to explain 
the excess emission at 4.5fim and 5.8/xm. AGNs and kpc- 
scale circumnuclear dust are commonly detected in nearby 
elliptical galaxies. In the Palomar spectroscopic survey of 
nearby galaxies, about half of the ellipticals show detectable 
emission-line nuclei, while most of which are classified as 
LINERs (Ho, et al. 1997). On the other hand, recent obser- 
vations from Hubble Space Telescope (HST) images show 
that circumnuclear dust appear in about ~ 40% of ellipti- 
cals in their optical images (Tran, et al. 2001; Lauer, et al. 
2005; Simoes Lopes, et al. 2007). PAHs were thought to be 
rare in ellipticals considering the sputtering destruction in 
the hot plasma environment. However, recently, Kaneda et 
al.(2008) reported detection of PAHs emission features in 
14 out of 18 dusty ellipticals, implying that PAHs are more 
common than ever thought in these systems. In the follow- 
ing discussion, we confine our interest to the role of different 
components playing in produce the excess infrared emission. 

Columns (4) and (5) of Table 1 summarize the classifi- 
cation on nuclear activity and circumnuclear dust morphol- 
ogy in optical images of our sample galaxies. All infrared 
core galaxies, both point-like and extended sources, show 
AGN activities, 8 of them are LINERs, except for NGC 4125, 
which is classified as a Transition object. Thus infrared core 
galaxies account for 41% of AGNs in our sample. With re- 
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spect to circumnuclear dust, 17 galaxies, about hall ol our 
sample, have been detected ol circumnuclear dust in optical 
band, all infrared core galaxies belong to this group. There- 
fore, the lact that both AGNs and circumnuclear dust coin- 
cide with a central infrared core makes it difficult to clarity 
their contributions to the infrared excess emission. 

Optical observations ol early-type galaxies show that 
active early-type galaxies more tend to possess circumnu- 
clear dust than inactive ones (Tran et al. 2001; Krajnovic 
& Jaffe, 2002; Xilouris & Papadakis, 2002; Simoes Lopes et 
al. 2007). Simoes Lopes et al. (2007) reported that 100% 
AGNs show ciucumnuclear dust leature in optical image, in 
contrast, only 2 out ol 15 ellipticals are detected with dust. 
This correlation can be understood in two ways. Firstly, the 
accretion ol a black hole requires luel, and circumnuclear 
dust is a good indication ol cold gas inflow and thus an ev- 
idence ol luel supply lor the central ol AGN. Nevertheless, 
the existence ol ellipticals without visible dust but hosting 
AGNs, although in small amount: NGC 2832, NGC 3608, 
NGC 5982, exclude the presumption that nuclear dust is 
necessary for nuclear activity. On the other hand, Temi et 
al. (2007a, b) suggested that AGN feedback may play an 
important role in transporting dust from central reservoir 
to the interstellar space in ellipticals, this mechanism may 
also assist to establish such a correlation. 

Figure 5 shows plots ol 60/j.m luminosity versus 100/im 
luminosity based on IRAS data (Knapp et al. 1989), most 
ol the lar-infrared detected galaxies are infrared core galax- 
ies hosting luminous AGN, while 4 inactive galaxies and 1 
non-core AGN also show detection but lower luminosities on 
the whole. This is consistent with the lact that AGN nuclei 
preferentially exist in optically dusty ellipticals. Although 
for the most powerful AGNs: NGC 1052, NGC 315, NGC 
4261, the lar-infrared emission may come from strong radio 
continuum powered by AGN, it is unlikely that this is the 
dominant contribution for other lainter AGNs, since their 
lar-infrared emission are too strong comparing with the ex- 
trapolation ol radio continuum. Additional evidence arises 
from resolved lar-infrared emission in NGC 5077 (Temi et 
al. 2007b) and equal dust mass between that inferred from 
lar-infrared measurement and that from optical extinction 
observation in NGC 4125 (Bregman et al. 1998). In dusty el- 
lipticals the dust mass inferred from lar-infrared is generally 
an order ol magnitude greater than that from optical extinc- 
tion observation (Goudfrooij et al. 1994), indicating that the 
distribution ol dust is extended. Therefore, the higher lar- 
infrared luminosities in luminous AGNs than laint AGN and 
inactive galaxies reveal that correlation between AGN and 
dust does not only hold for dust at circumnuclear region but 
also for extended dust throughout the galaxy. 

Figure 6 shows the 15 GHz nuclear radio luminosity ver- 
sus nuclear [O777]A5007 luminosity for all AGNs in our sam- 
ple. The nuclear 15 GHz radio luminosities are mainly taken 
from Nagar et al. (2005), and the nuclear [O777JA5007 lumi- 
nosities are taken from Ho et al.(1997), which have been cor- 
rected for reddening effect by using the Balmer decrement. It 
is obvious that infrared core AGNs distinctly separate from 
non-core AGNs with higher luminosities both at optical and 
radio wavelength. Especially, for [O777]A5007 emission, core 
AGNs and non-core AGNs are completely distinguished at 
log(Lo[iii]) ~ 38.5. With respect to radio emission, 7 out 
ol 9 core galaxies have been detected ol compact radio nu- 
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Table 3. Total and Excess Infrared Flux Densities. 



3 



Galaxy 


5T3.6 


■St4.5 


ST5.8 


Sts.q 


ST24 




-SE5.8 


■Se8.0 








Name 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


"3.6-8 


Q!8-24 


a radio 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 315 


207.7 ±20.7 


122.7 ±12.2 


74.0 ±7.4 


51.4 ±5.1 


100.1 ± 10.0 


1.3 ±0.3 


4.4 ±0.5 


10.0 ±0.6 


2.5 ±0.1 


2.1 


0.27 


NGC1052 


337.7 ± 33.7 


207.0 ± 20.7 


163.4 ± 16.3 


136.1 ± 13.6 




17.6 ±0.5 


38.7 ±0.5 


65.9 ±0.4 


1.8 ±0.1 


1.3" 


0.33 


NGC4261 


355.6 ±35.5 


209.1 ±20.9 


134.6 ± 1.3 


83.8 ± 8.3 


50.5 ± 5.0 


0.8 ±0.1 


2.0 ±0.2 


3.9 ± 0.3 


2.0 ±0.1 


2.3 


-0.24 


NGC5322 


360.2 ± 36.0 


216.2 ±21.6 


146.6 ± 14.6 


90.6 ±9.0 


33.5 ±3.3 


1.3 ±0.4 


4.0 ±0.6 


8.4 ±0.4 


2.3 ±0.5 


1.2 


0.15 


NGC3226 


117.1 ± 11.7 


61.7 ±6.1 


44.1 ±4.4 


33.5 ± 3.3 


30.1 ±3.0 


1.4 ±0.2 


4.6 ±0.3 


10.2 ±0.4 


2.3 ±0.2 


1.0 


0.11 


NGC4125 


541.2 ± 54.1 


309.2 ± 30.9 


211.0 ± 21.0 


122.5 ± 12.2 


31.9 ±3.2 






4.5 ±0.8 




1.8 




NGC4278 


386.3 ±38.6 


233.2 ±23.3 


149.7 ± 14.9 


110.6 ±11.0 


29.1 ± 2.9 






11.0 ±2.0 




0.9 


0.29 


NGC4374 


773.6 ± 77.3 


462.9 ± 46.3 


284.9 ±28.4 


180.4 ± 18.0 


28.9 ±2.8 






6.7 ± 1.3 




1.3 


0.13 


NGC5077 


162.6 ± 16.3 


94.1 ±9.4 


79.9 ±8.0 


37.8 ± 3.7 


21.5 ± 2.1 






1.9 ± 0.5 




2.2 





Notes — Col(2) — (6):Total infrared flux density at 3.6,4.5,5.8,8.0 and 2Afim. Col(7) — (9):Flux density of central excess emission at 3.6, 
5.8, 8.0/im. Col(10): Power-law index of central excess emission between 3.6 — 8/im, derived from 3.6, 5.8 and 8.0fim through method 

described in Section 3. Col(ll): Power-law index of central excess emission between 8 — 24/im, derived from Col(6) and Col(9). Col(12): 
Spectral index of nuclei derived from flux densities at 5GHz and 15GHz, which are taken from Nagar et al. (2005). 

a Since NGC1052 has no MIPS data, the power-law index is derived from the mid-infrared spectra, that is given by Kaneda et al, 2008. 



Table 4. Total and Excess Infrared Luminosities. 

4 



Galaxy £iaT3.6 L vTi5 £i/,T5.8 £^,T8.0 -£<iaT24.0 £iaE3.6 ^v,Eb.S £i/,E8.0 

Name (10 2s erg / sec) (10 28 erg / sec) (10 2s erg/sec) (10 28 erg/ sec) (10 28 erg/see) (10 28 erg/see) (10 28 erg/sec) (10 28 erg/see) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC315 


107.61 


63.57 


38.35 


26.64 


51.28 


0.68 


2.26 


5.08 


NGC1052 


12.80 


7.85 


6.19 


5.16 




0.67 


1.47 


2.50 


NGC4261 


52.41 


30.82 


19.84 


12.35 


7.41 


0.11 


0.35 


0.57 


NGC5322 


43.03 


25.83 


17.52 


10.82 


3.98 


0.16 


0.47 


1.00 


NGC3226 


7.67 


4.04 


2.89 


2.20 


1.95 


0.09 


0.30 


0.67 


NGC4125 


37.92 


21.67 


14.78 


8.58 


2.24 






0.31 


NGC4278 


4.35 


2.63 


1.69 


1.24 


0.33 






0.12 


NGC4374 


26.12 


15.63 


9.62 


6.09 


0.97 






0.22 


NGC5077 


32.07 


18.55 


15.71 


7.45 


4.27 






0.38 



Notes — Col(2) — (5):Total infrared luminosity at 3.6, 4.5, 5.8, & 8.0fim. Col(6) — (8):Luminosity of central excess emission at 3.6, 5.8, 

8.0/^m. 



clei, all of which exhibit flat spectrum — as shown in the last 
column of Table 3 — and high brightness temperature with 
T > 10 7 K indicating non-stellar origin. 

Figure 6 strongly suggests that the activity of cen- 
tral black hole accounts for infrared excess emission. If in- 
frared excess luminosity is correlated with AGN luminosity, 
a infrared-red core will intrinsically exist for every AGN, 
but be too weak to be detected in cases of low luminosity. 
The four panels of Figure 7 show 8.0pm excess luminos- 
ity plotted against nuclear 15GHz radio luminosity, MIPS 
24pm luminosity, nuclear [O///]A5007 luminosity and nu- 
clear X-ray luminosity at 2-10keV. The linear fit of data 
points yields linear-correlation R values of 0.65, 0.92, 0.32 
and 0.78, respectively. While the correlation between excess 
8pm emission and 24pm emission can be simply attributed 
to similar origins, correlations between 8pm excess and radio 
and hard X-ray emission, further support an AGN origin for 
infrared core. In comparison, the weak correlation between 
excess emission and [O77/]A5007 luminosity might be re- 
sulted from both extinction and shock ionization. A possible 
interpretation for the considerable scatter in Figure 7 is un- 
resolved non-AGN sources. Infrared core galaxies with weak 
excess emission in Figure 3 generally have extended struc- 
ture, while higher- luminosity sources mainly display point- 



like morphology, suggesting non-AGN contribution in low 
luminosity sources. Particularly, the weak emission at 8pm 
in extended sources imply PAHs origin, which are generally 
detected in dusty ellipticals (Bregman et al. 2006; Kaneda 
et al. 2008). Yet, PAHs emission, even if exist (For example: 
NGC 1052, Kaneda et al. 2008), will still play a minor role 
in strong and point-like sources considering the remarkable 
emission at wavelength shorter than 8pm in these objects. 

The contribution of PAHs can be further evaluated by 
comparing spectral indices at long and short wavelength. 
Since uncertainties of excess emission are serious at wave- 
length shorter than 5.8pm, here we only consider three 
galaxies (NGC 315, NGC 1052, NGC 3226) given their re- 
liable measurement at short wavelength. The approach to 
derive excess flux density described in Section 3 is based 
on the assumption that excess emission through IRAC band 
satisfy a single power-law distribution with a uniform alpha. 
In order to compare spectral indices between different IRAC 
bands, the flux density need to be reevaluated in another 
way. For NGC 315 and NGC 3226, the proportion of non- 
stellar component at 3.6pm flux density is negligible, thus 
we assume that all 3.6pm emission is due to old stellar pop- 
ulation, following this assumption the derived spectral in- 
dex a 4 . S -5.8 is 2.0 for NGC 315 and 3.4 for NGC 3226, 
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Figure 7. Comparison between 8/im excess emission and a set of AGN indicators. Open circles denote point-like sources and filled 
triangles denote extended sources. Y-axis refers to luminosities of excess emission at 8fim and x-axis refers to 15GHz core luminosities, 
MIPS 24/im luminosities, nuclear [O///]A5007 luminosities, and nuclear luminosities at 2-10keV, from left to right. The linear-correlation 
coefficient R is also shown. 
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Figure 6. Nuclear [O/iY]A5007 luminosity plotted against 
15GHz radio luminosity, symbols with red color refer to infrared 
core galaxies, stars are point-like core galaxies, triangles are ex- 
tended core galaxies, and crosses are other LLAGNs without in- 
frared core. 



while 05.8-8.0 = 2.8 for both galaxies. The spectral indices 
are similar enough to rule out the possibility of dominant 
PAHs emission. For NGC 1052, we simply derive the spec- 
tral index 03.6,4.5,5.8 following the same approach described 
in Section 3, the derived value is 1.7, also consistent with 

03.6,5.8,8.0 = 1-8. 

In normal elliptical galaxies, mid-infrared emission is 
dominated by circumstellar hot dust around AGB stars 



(Knapp et al. 1992; Temi et al. 2007a). Nevertheless, the 
tight correlation between MIPS 24pm emission and 8pm 
emission in Figure 7 suggests that the sources of both are 
the same in infrared core galaxies. To examine the origin of 
MIPS 24pm emission, we compare 24pm luminosities with 
optical B-band absolute magnitude Mb in Figure 8. As Temi 
et al. (2007a) noted in their work, for non-active ellipticals, 
Mb scales with 24pm luminosity, supporting a stellar origin 
for 24pm emissions. Half of galaxies with 8pm excess emis- 
sion also show 24pm excess with respect to normal galaxies, 
while the other four galaxies do not show obvious signs of 
24pm excess. Thus the tightness of the correlation between 
MIPS 24pm emission and 8pm emission should be taken 
with care when it comes to faint sources, which could be 
diluted by circumstellar emission. 

Three AGN-associated mechanisms can be responsible 
for the observerd excess infrared emission. The first is ther- 
mal emission from nuclear hot dust heated by central AGN, 
which could be expected in normal Seyfert galaxies. As dis- 
cussed earlier, in active ellipticals, optically observed cir- 
cunmuclear dust feature is common , generally with scales 
of a few hundreds pc. With such distance from nuclei, dust 
could not be heated sufficiently by AGN to generate ob- 
served emission at IRAC band (van Bemmel et al. 2003 & 
2004). Hot dust responsible for observed emission, if exist, 
will lie within a central region of several to about ten pc, 
where is traditionally assumed as torus and is difficult to 
be resovled. However, the viewpoint that LLAGNs possess 
torus as Serferts is challenged by the observed low X-ray col- 
umn densities in LLAGNs (Terashima et al. 2002; Gonzlez- 
Martn et al. 2006) and high detection rate of optical compact 
core (Chiaberge et al. 1999), which indicates mild obscura- 
tion. The disappearance of dust emission is also predicted 
by "wind-torus" model, in which obscuring torus is essen- 
tially clumpy dusty wind emanating from accretion disk and 
could not maintain while accretion rate declines to the level 
insufficient for out flow (Emmering et al. 1992; Konigl & 
Kartje, 1994; Elitzur et al. 2006). Without hot dust, a sec- 
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Figure 8. MIPS 24fim luminosity plotted against optical B-band 
absolute magnitude. Symbols have the same meanings as in Fig- 
ure 5 and Figure 6. 



ond choice to interpret the infrared excess is the synchrotron 
emission from thermal electrons in accretion flows, existing 
models of RIAF predict a submillimctcr to infrared bump in 
the SED (Quataert et al. 1999; Yuan et al. 2003). Finally, it 
is worthwhile to notice that 8 out of 9 core galaxies in our 
sample present compact radio core, and 4 of which (NGC 
315, NGC 1052, NGC 4261, NGC 4374) belong to the cat- 
egory of FR I radio galaxies with kpc scale jets. There is a 
long-time consideration that the SED of FR I radio galax- 
ies, or even other types of LLAGNs, could be dominated by 
synchrotron emission from the jet (Chiaberge et al. 1999; 
Yuan et al. 2002; Falcke et al. 2004). One typical example of 
jet dominated LLAGN is M87 (Shi et al. 2007; Perlman et 
al. 2007), where emission from bright knots of jet could be 
seen at all IRAC bands and nuclear mid-infrared emission is 
primarily due to the non-thermal emission from the base of 
the jet. Yet, unlike M87, none of core galaxies in our sample 
show any signs of jet emission at IRAC band, and our core 
galaxies generally have radio luminosities lower than M87 
by a factor of 100, but comparable infrared excess emission. 

As shown in Table 3, for 5 sources in our sample with 
detectable excess emission throughout the IRAC band, we 
give their spectral indices «3.6-8 following the method de- 
scribed in Section 3. The spectral indices of the 5 galaxies 
are similar, with values ~ 2, implying a common origin of 
infrared excess. Such an index value is too steep for a jet 
dominated SED, which is not expected to be significantly 
larger than 1.0 (Markoff et al. 2003; Perlman et al. 2007), 
and too flat for the model of radiation from RIAF, which 
lies around 3 (Yuan et al. 2003 and private communication 
with Yuan). A combination of this two components might 
interpret our derived indices, which could be expected in FR 
I radio galaxies (Wu et al. 2007). However, it should be no- 
ticed that the emission characters at radio wavelength of this 
5 core galaxies are not uniform, while 3 galaxies (NGC 315, 
NGC 1052, NGC 4261) are FR I galaxies with large scale jets 
and strong radio emission, NGC 3226 and NGC 5322 only 
show relatively faint compact core and show no extended ra- 
dio emission (Nagar et al. 2005). Thus it is questionable why 



the galaxies with distinctive radio-loudness could produce 
similar spectral energy distribution through hybrid emis- 
sion including jet and accretion flows. On the contrary, this 
value is consistent with that of general spectral index of 
Seyfert nuclei (Alonso-Herrero et al. 2003), supporting ori- 
gin of thermal emission from hot dust. In addition, Barth et 
al. (1999) detected polarized broad H a emission in and only 
in the 3 FR I galaxies mentioned above (NGC 315, NGC 
1052, NGC 4261), out of 14 LLAGNs, and all of them also 
show considerable absorption X-ray columns (~ 10 22 cm' 2 ) 
(NGC 1052: Guainazzi et al. 2000; NGC 4261: Sambruna 
et al. 2003; NGC 315: Gonzlez-Martn et al. 2006 ), indicat- 
ing the existence of obscuring structure and at least some 
contribution of thermal emission. Since these three objects 
have relatively higher AGN luminosities, it is possible they 
more closely resemble classical AGNs than other objects in 
our sample. Therefore, while we are more inclined to at- 
tribute the infrared excess of relatively higher luminosity 
AGNs to dust emission, for evidence of torus in some of 
them, and for their similar spectral indices with Seyfert2s, 
the infrared emission mechanism of fainter objects is more 
uncertain. Our results show that infrared excess emission 
decreases with the decrease of other AGN indicators, and 
current data show no sign for a change of infrared emission 
mechanism. On the other hand, although LLAGNs generally 
have low X-ray column density, X-ray absorption does not 
have a direct connection with infrared emission. Firstly, the 
X-ray inferred column density is determined by obscuring 
material along our line of sight, while the infrared emission 
is integrated emission emission from all sources surrounding 
AGN. Secondly, while both produce X-ray absorption, the 
contribution of dust-free gas is larger as compared to dusty 
gas. The important role of dust-free gas in X-ray absorption 
is supported by lower column density inferred from redden- 
ing effect than that from X-ray absorption. (Maccacaro et 
al. 1982; Maiolino et al. 2001). The insignificant X-ray ab- 
sorption does not automatically indicate disappearance of 
dust emission. By assuming Galactic dust-to-gas ratio and 
"standard" model of interstellar dust, the column density 
N H ~ 2 x 10 21 t v cm" 2 (Elitzur 2008), a neutral hydrogen 
column density with a few 10 21 cm~ 2 could still cause a con- 
siderable absorption in optical/UV band and thus thermal 
re-emission in infrared band. 

After all, we should acknowledge that current data are 
insufficient to draw a ultimate conclusion. At the resolution 
of IRAC, it is difficult to further identify the mechanism to 
produce nuclear infrared emission among different possibil- 
ities. Comprehensive study of the SED might provide more 
information on this issue. 



5 CONCLUSION 

We performed the Spitzer IRAC observations of 36 local el- 
liptical galaxies. 9 out of 36 galaxies display red core struc- 
ture with nuclear infrared excess emission. The infrared ex- 
cess emissions only and universally appear in galaxies with 
relatively luminous central AGN, strongly support a relation 
between the two. We also confirmed the correlation between 
the activities of AGN and optically observed circumnuclear 
dust and found this correlation also holds for extended dust 
in elliptical galaxies. We found correlation with considerable 
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scatter between the luminosity of central AGN and excess 
emission, which indicates unresolved non-AGN contamina- 
tion of excess emission in low-luminosity sources. While the 
specific mechanism to produce infrared emission could not 
be identified by current data, thermal origin from hot dust is 
supported by similar infrared spectral indices with Seyfert 
galaxies. In order to clarify origins for LLAGNs' infrared 
emission, a further study calls for multi-wavelength study, 
which will be provided in our next work. 
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